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Abstract 

We experimentally investigate the nucleation of surface nanobubbles on PFDTS-coated silicon as 
a function of the specific gas dissolved in the water. In each case we restrict ourselves to equilibrium 
conditions (c = 100%, Tuq^id — Tsubstrate)- Not only is nanobubble nucleation a strong function of 
gas type, but there also exists an optimal system temperature of ^ 35 — 40 where nucleation is 
maximized, which is weakly dependent on gas type. We also find that contact angle is a function 
of nanobubble radius of curvature for all gas types investigated. Fitting this data allows us to 
describe a line tension which is dependent on the type of gas, indicating that the nanobubbles are 
sat on top of adsorbed gas molecules. The average line tension was r ^ — 0.8nN. 
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I. INTRODUCTION 



Surface nanobubbles are nanoscopic gaseous domains that exist on the sohd/hquid in- 
terface p]-[TT]. They are surprisingly stable [TOl [T214T4] . surviving orders of magnitude 
longer than the classical diffusive life time. Since their recent discovery, several key observa- 
tions have been made regarding the 'ideal' conditions for nucleation p^SHTTj . These include 
the (almost p!8l [19j ) uniqueness for the formation in water, or in solutions containing at 
least some part water. On the other hand, changing the substrate chemistry can lead to 
differences in the typical nanobubble sizes and distributions [6l [20l-[22j. but nanobubbles 
are always found as long as sufficient care is used as well as an appropriate experimental 
technique/pretreatment [SHU [6l [23l - [26] . 

One area that has seemingly escaped detailed research is the dependence of nanobub- 
ble formation on the type of gas that is dissolved in the water. Almost all studies to date 
investigate nanobubbles made from 'air', with the few exceptions including carbon diox- 
ide [121 HQ EZ], butane [26l [28], and nitrogen [26j (we omit here the hydrogen and oxygen 
nanobubbles created using electrolysis [29l - [3Tl| . since gas was not necessarily dissolved in 
the water). Of these specific gases, carbon dioxide and butane seem peculiar choices. Both 
of these gases have triple points above room temperature so can liquefy if the pressure is 
sufficiently high, i.e. if they fill a nanobubble which is sufficiently small. As an example, 
carbon dioxide was used in Reference [12j because of its strong rotational fine structure 
signal, which allowed use of the attenuated total reflection infrared spectroscopy technique 
to show conclusively that nanobubbles contain gas. However the critical pressure and tem- 
perature of CO2 are Pc = 7.4 MPa and Tc = 31.0 °C |32j, respectively, so we would expect 
liquefaction to occur in these C02-filled nanobubbles if their radii of curvature was lower 
than Rc = 2^/{pc — Po) ^ 20 nm (note that the nanobubbles in Reference [12j were much 
larger than this). In butane-saturated water, nanodroplets have already been reported on 
Si(lOO) by Kameda et al [21j. 

It is clear that gas type plays an important role, but no study has ever been carried 
out to investigate the effects of different gases. Thus, we have directly investigated the 
role of gas type on the formation of surface nanobubbles. We chose a selection of different 
gases, including noble gases (helium and argon), diatomic gases (hydrogen, nitrogen, and 
oxygen), and more complex gases (carbon dioxide and methane). Furthermore, for each gas 
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type we restricted ourselves to systems where the hquid was saturated with the gas and 
specificaUy not supersaturated, foUowing recent observations that supersaturation is not a 
requirement for nanobubble formation [17j. Thus our observations provide information on 
surface nanobubble nucleation in equilibrium. 

II. EXPERIMENTAL METHOD 

The substrate was a silicon wafer that had been hydrophobised with a self assem- 
bled monolayer of perfluorodecyltrichlorosilane (PFDTS), following the guidelines of Ref- 
erence [331134]. In brief, we degassed the chamber that contained the uncoated substrate to 
a pressure of ^ 50 /xbar (i.e. below the vapor pressure of the PFDTS). Next, we opened this 
chamber to a reservoir of degassed PFDTS for 5 min so that the silane molecules adsorbed 
onto the substrate. Finally, we closed the system off to the PFDTS reservoir and opened 
it to a reservoir of degassed pure water for half an hour, hence increasing the pressure of 
the system to the vapor pressure of water to allow the reagents to react. The resulting rms 
roughness was 0.4 nm, and the contact angles were ^ 110^ (equilibrium), ^ 116° (advanc- 
ing), and ^ 104° (receding). This surface has no known phase transitions in the temperature 
range investigated in the present study. 

The substrate was mounted on a temperature controlled sample plate (331 temperature 
controller, Lakeshore, USA) before a purpose-built atomic force microscope (AFM) liquid 
cell was firmly mounted on top. The AFM was an Agilent 5100 operated in tapping mode 
with a typical scan speed of 4 — 5 /xm/s. The AFM cantilevers were hydrophilic, Au-back- 
coated Si3N4 Veeco NPG probes, with typical spring constants and resonance frequencies in 
water of 0.58 N/m and 25 - 35 kHz. We operated the AFM at a set point of 90%. 

For the liquid we used ultrapure water (Simplicity 185 purification system, Millipore SAS, 
France). This was placed within a stainless steel container that included a pump-out port, 
a re-pressurization port, and a digital pressure gauge. 

Our study was to investigate the role of the gas type on nanobubble nucleation, so we now 
explain the procedure that we used to prepare the gas dissolved in the water. The first stage 
of preparation was to thoroughly degas the water to remove as much air as possible. During 
this stage the pressure was reduced to < 20mbar and, once this pressure was reached, we 
continued to pump for at least a further ^ 30mins. The container of water was then closed 
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off to the pump and placed within a heat bath which was at the correct temperature for the 
experiment at hand. 

The stainless steel container had a volume of 100 mL but we only treated 40 mL of water 
at a time. This allowed us to repressurise the evacuated volume above the water level with 
a specific type of gas. The number of molecules of a specific type that dissolve in water 
at equilibrium is well known [32j, thus the magnitude of the repressurisation was to latm 
plus the required number of molecules for dissolution multiplied by kT/V^^ where kT is the 
thermal energy and is the 40 mL volume of water. We plot the solubilities of the different 
gases used in this study, as functions of temperature, in [T] After the container had been 
pressurized we closed it off to the gas supply and maintained the system's seal for 2 — 3 
gas-diffusion time scales. 

Once the liquid was prepared it was injected into the liquid cell of our atomic force 
microscope using a syringe pump. The syringe was at the same temperature as both the 
water and the substrate. The AFM was placed inside a glass environment-control chamber 
that had had its air displaced with one atmosphere of the same gas type that was dissolved in 
the liquid. Thus, any gas exchange between the liquid and the environment would not lead 
to cross contamination of the gas type. This means that the gas content of the water was 
maintained throughout the measurements and our results provide information of nanobubble 
nucleation in fully equilibrated systems. 

III. RESULTS AND DISCUSSION 

We present typical images of nanobubbles nucleated from different gases at 25 in [2} 
Each image represents 2 /xm x 2 /xm, but the height scales differ as described in the caption. 
Each gas type resulted in distinct nanobubble sizes/densities, with methane being the only 
gas that led to zero nucleation at this temperature. Note that each set of experimental 
conditions was repeated in their entirety, including the water preparation stage, and for a 
given experiment several areas of the substrate were scanned. 

To quantify the differences between the gas types we plot the average radii of curvature 
and diameters of nanobubbles as functions of gas type in [2)i and [2^, respectively, as well 
as the calculated total volume of nanobubbles on a 2 jam x 2 /xm area in [3} where we plot 
the 'gas type' along the abscissa in order of increasing molecular weight and solubility (note 
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FIG. 1: Concentration of saturated gases in water as a function of temperature. Colors are for H2 
(dark green), He (blue), CH4 (red), N2 (green), O2 (magenta), Ar (cyan), and CO2 (black). Note 
that the data for H2 and He overlie each other. 



that each data point is an average of several different areas and several experimental runs 
- we are displaying the average value as a 'standard' measure of nanobubble populations, 
see References pTi |35j). If solubility was the control parameter governing nanobubble nu- 
cleation, we would expect the data in |3] to be a monotonically varying function. Clearly 
this is not the case, with no obvious functional form able to describe the data. This means 
that solubility is not the governing factor for nanobubble nucleation in a system that is 
in equilibrium (i.e. with the substrate and liquid temperatures being equal, and the gas 
saturated within the liquid). 

Secondly, although adsorption strengths of the various gases are not known for our sub- 
strate, we expect the ordering of these to be He and H2 (weak); Ar, O2, and N2 (medium); 
CH4 and CO2 (strong) [361438]. Clearly there is no direct dependence on this either. Hence 
we can rule out the possibility of nanobubbles solely nucleating from the bulk desorption of 
dense adsorbates, which may exist in the form of micropancakes [34j. 

Another possible origin of the nucleation may be the excess gas that swells the density- 
depleted layer immediately next to an immersed hydrophobic solid [39l HOj. However, the 
dependency of the amount of swelling on gas type is currently disputed [391 SD 112] • Thus 
none of the three most likely possible causes for nanobubble nucleation, namely formation 
from gas dissolved in the bulk, formation from molecules adsorbed on the substrate, or 
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FIG. 2: Typical images of nanobubbles formed from different gas types at 25 ^C. Gases are (a) 
H2, (b) He, (c) CH4, (d) N2, (e) O2, (f) Ar, and (g) CO2. Images are 2 /im x 2 /xm; height scales 
are (b,d) 20 nm and (a,c,e-g) 250 nm. (h) Average radii of curvature and (i) average diameters of 
nanobubbles versus gas type (error bars correspond to one standard deviation). 

formation from the gas-enriched layer near the hydrophobic substrate, can be individually 
responsible. 

In order to try to understand the gas dependency of nanobubble nucleation further, we 
proceeded to investigate the effects of different temperatures on nucleation. For this, we re- 
peated the experiments for methane, nitrogen, and oxygen, this time at system temperatures 
of 30 ^C, 35 ^C, 40 ^C, and 45 ^C. 

Typical images of nanobubbles created from the three g ases, as functions of system tem- 
perature, are shown in |4} It is evident that an optimal temperature exists for each type of gas 
that leads to maximum nanobubble nucleation. As an example, no nanobubbles nucleated 
from the methane-saturated water at 25 ^C, but nucleation occurred for every temperature 
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FIG. 3: The total nanobubble volume on a 2 x 2 area as a function of gas type at a system 
temperature of 25*^0. The abscissa is ordered with increasing molecular weight and solubility. If 
gas solubility was the control parameter governing nucleation we would expect the data in this 
graph to vary monotonically. 

we tested above this, with maximum production at 40 ^C. These optimal temperatures are 
more clearly visible in the plots of the total nanobubble volume in [5j Other work pTl l35] 
has described a similar maximum in azr-fiUed nanobubbles as being due to a maximum in 
solubility. However, if we compare the curves of [5] to the temperature dependencies of the 
solubilities in [T} it is clear that there is no corresponding peak in solubility, i.e. we reiterate 
that solubility is not the control parameter for nanobubble nucleation in gas-equilibrated 
systems. 

In total, over 300 nanobubbles were imaged for the present study, so we have sufficient 
statistics to extract information for the contact angle versus radius of curvature. We plot 
this for methane, nitrogen, and oxygen in [6| where we use different symbols and colors to 
represent the differing system temperatures. For all of the gases in the present study the 
contact angle was found to increase with increasing radius of curvature, before leveling off 
to^^- 180^asi?^oc. 

When contact angle is a function of the radius of curvature, it is necessary to introduce 
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FIG. 4: Temperature dependence of nanobubbles formed from methane, nitrogen, and oxygen- 
saturated water at temperatures of (a) 25 ^C, (b) 30 ^C, (c) 35 ^C, (d) 40 ^C, and (e) 45 ^C. Images 
are 2 /xm x 2 /xm. Height scales are (methane, 40*^0) 230 nm, (methane, other) 40 nm; (nitrogen) 
40 nm; (oxygen, 25 ^C) 250 nm, (oxygen, other) 40 nm. 
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FIG. 5: The total nanobubble volume on a 2 jam x 2 jam area as a function of gas type and 
temperature for methane (triangles), nitrogen (circles), and oxygen (asterisks). There is a clear 
maximum in density between 35 and 40 for all three gas types. 
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a line tension. Brenner and Lohse [lOj suggest a functional form for this effect of 

cos 6>oo - cos 6>o 

cos (9 = cos (900 , , j^,. . (1) 

1 + n/o 

Here, ^oo is the value of the contact angle at large radii, 6^ is at diminishingly small radii, 
and 5 is the length scale over which we expect line tension to take effect. Hence, we now 
choose to use the functional form of Reference for the nanoscale correction to contact 
angle demonstrated in |6| For the three fitting parameters, ^oo, ^o, and 5^ we set the two 
angles to 180^ and 90^, respectively, and use 5 as the single fitting parameter for the data. 
These choices for the angles were selected because (i) it is clear from the data in |6]that the 
contact angle tends to 180^ with increasing radius of curvature, and (ii) the contact angle 
rapidly reduces for smaller bubbles but we do not want or expect the nature of the material 
to alter from hydrophobic to hydrophilic with decreasing radius of curvature. Within this 
framework line tension would be r = — 7(^, where 7 is the surface tension of water. A typical 
fit is shown in[7| where we present the 35 data for oxygen. In this case, a line tension of 
r — 0.6nN becomes important at a length scale oi 5 9nm. 

Values of 5 for each of the different gases used in this study, as well as for the different 
temperatures used for methane, nitrogen, and oxygen, are presented in |8| There is a 
large spread in 5 with gas type, which gives insight into the possible configuration of the 
nanobubble-substrate geometry. Different values of 8 for different gas types means that 
the contact angle is different for different gas types. This is puzzling since contact angle 
is dependent on the three differing surface energies of the three different interfaces (solid, 
liquid, and gas), but it is always the densest phase that contributes the most to these values. 
(As an example, degassing the air above the water level in a glass beaker does not lead to 
a change in contact angle of the meniscus - the air has been evacuated, but the water and 
glass remain.) Thus, the specific type of gas should not have a noticeable effect on contact 
angle because this is dominated by the water and solid. The way to resolve this issue of 
variable contact angle with gas type would be if the nanobubble was sat on top of a dense 
adsorbate ("micropancake," [43j) of gas molecules. As pointed out by Reference [llj, this 
is the most probable configuration for a nanobubble - if a nanobubble is sat on top of a 
dense adsorbate, it is the binding energy of gaseous molecules to adsorbed molecules that 
should be considered and not the binding energy to the underlying solid. Not only does this 
explain the essentially constant value of contact angle for air- filled nanobubbles, regardless 
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FIG. 6: Contact angle as a function of nanobubble radius of curvature for (a) methane, (b) nitrogen, 
and (c) oxygen nanobubbles at 25 (circles), 30 "^C (right triangles), 35 (diamonds), 40 
(squares), and 45*^0 (up triangles). Dependence of contact angle on radius of curvature indicates 
that line tension must be considered. Data has been corrected for cantilever tip distortion. 
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FIG. 7: Fit of the 35 data for contact angle versus radius of curvature for oxygen. Fit is of the 
form cos9 = cos6>oo - (cos6>oo - cos6>o)/(l + R/S), with O^c = 0^, 9o = 90^, and 5 = 9nm. The 
corresponding hne tension is r = — 7^ = — 0.6nN. 

of the substrate, but it also explains the difference of contact angle with gas type. Of course, 
an underlying adsorbate would also provide a gas bank for the gaseous influx to balance the 
diffusive outflux in the dynamic equilibrium model of nanobubble stability [lOj. 

The mean value of line tension here, which we calculate as the average of the data in 
[8^,b, is r ^ — 0.8nN. Measurement of a negative line tension for nanobubbles/nanodroplets 
is expected[44j. Values in the literature include -3nN [27j, - -2.3 pN [45j, -0.2 nN [21j, 
and ^ -IpN [i44]. 

Returning to [8| we can again reiterate that we flnd no clear functional form linking 
solubility or adsorption strength to the different values of 6. Thus, we posit the following 
possible route to nucleation, using the noble gases. Argon should adsorb more strongly to 
the substrate than helium, resulting in less argon nanobubbles, but the opposite is found. 
The fact that there is approximately an order of magnitude more argon available in the 
bulk, however, should increase the density of the adsorbate, possibly to include multiple 
layers. Thus nanobubble nucleation may occur as a result of bulk desorption from the more 
weakly bound upper layers. For the diatomic molecules (hydrogen, nitrogen, and oxygen) 
it is a lot less clear. The strength of these adsorbates is dependent on the mean orientation 
of the molecules, and we would expect bulk desorption to occur much more readily than for 
the noble gases. The same is also true for the more complex molecules (carbon dioxide and 
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FIG. 8: Fitting parameter 6 as (a) a function of gas type at a temperature of 25 °C and (b) 
as a function of temperature for methane, nitrogen, and oxygen. The average hne tension is 
r -j6 - -0.8 nN. 

methane) . 

IV. CONCLUSIONS 

We have shown that gas type is a key parameter for the nucleation of nanobubbles. Not 
only do specific gases lead to the formation of more nanobubbles than others, but there also 
exists an optimal temperature for nanobubble nucleation between ^ 35 and ^ 40 ^C, 
which appears to be weakly dependent on gas type. 

Surprisingly, nanobubble nucleation is not directly dependent on either the solubility of 
the specific gas in water, or on the relative adsorption strength of the gas to the substrate. 
This indicates that nanobubbles do not form solely due to the amount of gas available in the 
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bulk, or from dense adsorbates (micropancakes) on the substrate. Furthermore, we should 
expect either zero or very small dependence on the specific gas type for the gas-enrichment 
layer thickness near a hydrophobic substrate, so this does not provide an adequate solution 
to nucleation either. 

Hence, nanobubble nucleation must come from a combination of several competing fac- 
tors. Certainly (at least) a two-stage nucleation process is most likely. On the one hand, 
more nanobubbles are found with increasing temperature for low temperatures, as would be 
consistent for systems reliant on desorption, whilst less nanobubbles are found with increas- 
ing temperature for high temperatures, now consistent with systems reliant on solubility. 
The cross over in behaviour is at a temperature of ^ 35 — 40 ^C. 

For all the nanobubbles thatwe investigated we found a dependence of contact angle 
on the radius of curvature. We introduced a line-tension term of the form proposed by 
Brenner and Lohse j^. On our PFDTS-coated silicon substrate, the average line tension 
was negative and equal to r ^ — 0.8nN. 
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